Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance (preprint)

The capacity to undertake whole genome sequencing (WGS) in public health laboratories (PHLs) has grown rapidly in response to COVID-19, and SARS-CoV-2 genomic data has been invaluable for managing the pandemic. The public health response has been further supported by the rapid upgrade and implementation of laboratory and bioinformatic resources. However, there remains a high degree of variability in methods and capabilities between laboratories. In addition to evolving methodology and improved understanding of SARS-CoV-2, public health laboratories have become strained during surges in case numbers, adding to the difficulty of ensuring the highest data accuracy. Here, we formed a national working group comprised of laboratory scientists and bioinformaticians from Australia and New Zealand to improve data concordance across PHLs. Through investigating discordant sequence data from Australia's first external SARS-CoV-2 WGS proficiency testing program (PTP), we show that most discrepancies in genome assessment arose from intrahost variation. While others could be remedied using reasonable, parsimonious bioinformatic quality control. Furthermore, we demonstrate how multidisciplinary national working groups can inform guidelines in real time for bioinformatic quality acceptance criteria. Provision of technical feedback allows laboratory improvement during a pandemic in real time, enhancing public health responses.

State-wide Genomic Epidemiology Investigations of COVID-19 Infections in Healthcare Workers: Insights for Future Pandemic Preparedness (preprint)

BackgroundCOVID-19 has resulted in many infections in healthcare workers (HCWs) globally. We performed state-wide SARS-CoV-2 genomic epidemiological investigations to identify HCW transmission dynamics and provide recommendations to optimise healthcare system preparedness for future outbreaks. MethodsGenome sequencing was attempted on all COVID-19 cases in Victoria, Australia. We combined genomic and epidemiologic data to investigate the source of HCW infections across multiple healthcare facilities (HCFs) in the state. Phylogenetic analysis and fine-scale hierarchical clustering were performed for the entire Victorian dataset including community and healthcare cases. Facilities provided standardised epidemiological data and putative transmission links. FindingsBetween March and October 2020, approximately 1,240 HCW COVID-19 infection cases were identified; 765 are included here. Genomic sequencing was successful for 612 (80%) cases. Thirty-six investigations were undertaken across 12 HCFs. Genomic analysis revealed that multiple introductions of COVID-19 into facilities (31/36) were more common than single introductions (5/36). Major contributors to HCW acquisitions included mobility of staff and patients between wards and facilities, and characteristics and behaviours of individual patients including super-spreading events. Key limitations at the HCF level were identified. InterpretationGenomic epidemiological analyses enhanced understanding of HCW infections, revealing unsuspected clusters and transmission networks. Combined analysis of all HCWs and patients in a HCF should be conducted, supported by high rates of sequencing coverage for all cases in the population. Established systems for integrated genomic epidemiological investigations in healthcare settings will improve HCW safety in future pandemics. FundingThe Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

Immune responses in COVID-19 respiratory tract and blood reveal mechanisms of disease severity (preprint)

Although the respiratory tract is the primary site of SARS-CoV-2 infection and the ensuing immunopathology, respiratory immune responses are understudied and urgently needed to understand mechanisms underlying COVID-19 disease pathogenesis. We collected paired longitudinal blood and respiratory tract samples (endotracheal aspirate, sputum or pleural fluid) from hospitalized COVID-19 patients and non-COVID-19 controls. Cellular, humoral and cytokine responses were analysed and correlated with clinical data. SARS-CoV-2-specific IgM, IgG and IgA antibodies were detected using ELISA and multiplex assay in both the respiratory tract and blood of COVID-19 patients, although a higher receptor binding domain (RBD)-specific IgM and IgG seroconversion level was found in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples was detected only when high levels of RBD-specific antibodies were present. Strikingly, cytokine/chemokine levels and profiles greatly differed between respiratory samples and plasma, indicating that inflammation needs to be assessed in respiratory specimens for the accurate assessment of SARS-CoV-2 immunopathology. Diverse immune cell subsets were detected in respiratory samples, albeit dominated by neutrophils. Importantly, we also showed that dexamethasone and/or remdesivir treatment did not affect humoral responses in blood of COVID-19 patients. Overall, our study unveils stark differences in innate and adaptive immune responses between respiratory samples and blood and provides important insights into effect of drug therapy on immune responses in COVID-19 patients.

Immune responses in COVID-19 respiratory tract and blood reveal mechanisms of disease severity (preprint)

Although the respiratory tract is the primary site of SARS-CoV-2 infection and the ensuing immunopathology, respiratory immune responses are understudied and urgently needed to understand mechanisms underlying COVID-19 disease pathogenesis. We collected paired longitudinal blood and respiratory tract samples (endotracheal aspirate, sputum or pleural fluid) from hospitalized COVID-19 patients and non-COVID-19 controls. Cellular, humoral and cytokine responses were analysed and correlated with clinical data. SARS-CoV-2-specific IgM, IgG and IgA antibodies were detected using ELISA and multiplex assay in both the respiratory tract and blood of COVID-19 patients, although a higher receptor binding domain (RBD)-specific IgM and IgG seroconversion level was found in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples was detected only when high levels of RBD-specific antibodies were present. Strikingly, cytokine/chemokine levels and profiles greatly differed between respiratory samples and plasma, indicating that inflammation needs to be assessed in respiratory specimens for the accurate assessment of SARS-CoV-2 immunopathology. Diverse immune cell subsets were detected in respiratory samples, albeit dominated by neutrophils. Importantly, we also showed that dexamethasone and/or remdesivir treatment did not affect humoral responses in blood of COVID-19 patients. Overall, our study unveils stark differences in innate and adaptive immune responses between respiratory samples and blood and provides important insights into effect of drug therapy on immune responses in COVID-19 patients.

Genomics-Informed Responses in the Elimination of COVID-19 in Australia (preprint)

BACKGROUND: A cornerstone of Australia’s ability to control COVID-19 has been effective border control, using an extensive supervised quarantine program. However, a rapid recrudescence in COVID-19 cases was observed in the state of Victoria in June 2020. Here, we describe the genomic findings that located the source of this second wave as a breach in supervised hotel quarantine and demonstrate the successful elimination of COVID-19 for a second time in Australia.METHODS: Genome sequencing was performed on all available SARS-CoV-2-positive samples in Victoria and integrated genomic and epidemiological investigation undertaken.RESULTS: At 31st January 2021, 20,451 COVID-19 cases were reported in Victoria; samples were sequenced from 75% of cases (15,431/20,451). Genomics revealed 98% (10,426/10,646) of locally-acquired cases during the second wave were derived from a single incursion from hotel quarantine, with the outbreak strain rapidly detected in other Australian states and territories. Phylodynamic analyses indicated an epidemic growth rate comparable to emerging variants, such as B.1.1.7 in the United Kingdom. Strict public health interventions resulted in the elimination of the outbreak strain by 29th October 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread.CONCLUSIONS: Rapid escalation of clonal outbreaks can occur from even a single breach of control practices, as revealed through our genomic ‘enhanced outbreak-detection' system. The subsequent elimination and rapid control of new SARS-CoV-2 incursions reinforce that decisive public health responses to emergent cases are effective even with high epidemic growth rates, and “elimination” should be favored in settings where this is achievable.FUNDING STATEMENT: The Microbiological Diagnostic Unit Public Health Laboratory (MDU PHL) and the Victorian Infectious Diseases Reference Laboratory (VIDRL) at The Doherty Institute are funded by the Victorian Government. This work was supported by the National Health and Medical Research Council, Australia (NHMRC); Partnership Grant (APP1149991), Investigator Grant to BPH (APP1196103), Investigator Grant to DAW (APP1174555), Research Fellowship to TPS (APP1105525), MRFF COVID-19 Genomics Grant (MRF9200006).DECLARATION OF INTERESTS: None to declare. ETHICS APPROVAL STATEMENT: Data were collected in accordance with the Victorian Public Health and Wellbeing Act 2008. Ethical approval was received from the University of Melbourne Human Research Ethics Committee (study number 1954615.3).

Tracking the COVID-19 pandemic in Australia using genomics (preprint)

BACKGROUND: Whole-genome sequencing of pathogens can improve resolution of outbreak clusters and define possible transmission networks. We applied high-throughput genome sequencing of SARS-CoV-2 to 75% of cases in the State of Victoria (population 6.24 million) in Australia. METHODS: Cases of SARS-CoV-2 infection were detected through active case finding and contact tracing. A dedicated SARS-CoV-2 multidisciplinary genomic response team was formed to enable rapid integration of epidemiological and genomic data. Phylodynamic analysis was performed to assess the putative impact of social restrictions. RESULTS: Between 25 January and 14 April 2020, 1,333 COVID-19 cases were reported in Victoria, with a peak in late March. After applying internal quality control parameters, 903 samples were included in genomic analyses. Sequenced samples from Australia were representative of the global diversity of SARS-CoV-2, consistent with epidemiological findings of multiple importations and limited onward transmission. In total, 76 distinct genomic clusters were identified; these included large clusters associated with social venues, healthcare facilities and cruise ships. Sequencing of sequential samples from 98 patients revealed minimal intra-patient SARS-CoV-2 genomic diversity. Phylodynamic modelling indicated a significant reduction in the effective viral reproductive number (Re) from 1.63 to 0.48 after the implementation of travel restrictions and population-level physical distancing. CONCLUSIONS: Our data provide a comprehensive framework for the use of SARS-CoV-2 genomics in public health responses. The application of genomics to rapidly identify SARS-CoV-2 transmission chains will become critically important as social restrictions ease globally. Public health responses to emergent cases must be swift, highly focused and effective.

Validation of a Single-step, Single-tube Reverse Transcription-Loop-Mediated Isothermal Amplification Assay for Rapid Detection of SARS-CoV-2 RNA (preprint)

2.IntroductionThe SARS-CoV-2 pandemic of 2020 has resulted in unparalleled requirements for RNA extraction kits and enzymes required for virus detection, leading to global shortages. This has necessitated the exploration of alternative diagnostic options to alleviate supply chain issues. AimTo establish and validate a reverse transcription loop-mediated isothermal amplification (RT-LAMP) assay for the detection of SARS-CoV-2 from nasopharyngeal swabs. MethodologyWe used a commercial RT-LAMP mastermix from OptiGene Ltd in combination with a primer set designed to detect the CDC N1 region of the SARS-CoV-2 nucleocapsid (N) gene. A single-tube, single-step fluorescence assay was implemented whereby as little as 1 L of universal transport medium (UTM) directly from a nasopharyngeal swab could be used as template, bypassing the requirement for RNA purification. Amplification and detection could be conducted in any thermocycler capable of holding 65{degrees}C for 30 minutes and measure fluorescence in the FAM channel at one-minute intervals. ResultsAssay evaluation by assessment of 157 clinical specimens previously screened by E-gene RT-qPCR revealed assay sensitivity and specificity of 87% and 100%, respectively. Results were fast, with an average time-to-positive (Tp) for 93 clinical samples of 14 minutes (SD {+/-}7 minutes). Using dilutions of SARS-CoV-2 virus spiked into UTM, we also evaluated assay performance against FDA guidelines for implementation of emergency-use diagnostics and established a limit-of-detection of 54 Tissue Culture Infectious Dose 50 per ml (TCID50 mL-1), with satisfactory assay sensitivity and specificity. A comparison of 20 clinical specimens between four laboratories showed excellent interlaboratory concordance; performing equally well on three different, commonly used thermocyclers, pointing to the robustness of the assay. ConclusionWith a simplified workflow, N1-STOP-LAMP is a powerful, scalable option for specific and rapid detection of SARS-CoV-2 and an additional resource in the diagnostic armamentarium against COVID-19. 3. Data summaryThe authors confirm all supporting data, code and protocols have been provided within the article or through supplementary data files.